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We study the dependence on filling and pressure of the superconducting and ferromagnetic 
critical temperatures of the ruthenocuprates, within the two-band model. At zero pressure, we find 
separate regions of coexistence of superconductivity and ferromagnetism as a function of filling, 
with contiguous regions merging together as pressure increases. As a function of pressure, a stronger 
enhancement of the magnetic phase results in a reduced pressure effect on the superconducting 
critical temperature. Comparison with recent experiments on the determination of the critical 
temperatures as a function of the pressure is also discussed. 
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I. INTRODUCTION 

The strong electron correlations in transition metal 
oxides are known to generate complex phase diagrams. 
Nonetheless, it was surprising to discover that one such 
oxide, the ruthenocuprate RuSr 2 GdCu20s (Ru-1212), 
shows a low-temperature phase in which homogeneous 
ferromagnetism and superconductivity coexist^* 2 . This 
new class of materials then adds itself to the metal- 
lic ferromagnets UGe2 (Ref. 0), ZrZn 2 (Ref. 0), and 
URhGe (Ref. H, all showing coexistence of unconven- 
tional superconductivity and itinerant ferromagnetism at 
low temperature (T ~ 1 K) and relatively high pressure 
(P-lGPa). 

The problem of the coexistence of ferromagnetism 
(FM) and superconductivity (SC) is puzzling due to the 
fact that there are at least two factors that would destroy 
superconductivity in a ferromagnetic medium: first, the 
exchange splitting lifts the energy degeneracy of the part- 
ners of a spin-up and spin-down Cooper pair; second, 
magnon exchange leads to repulsion for a singlet pair. 
Thus, the coexistence of ferromagnetic order and super- 
conductivity in the Ru-1212 and Ru-1222 compounds 
raises the question how these two antagonist states of 
matter can accommodate each other. Do both states 
coexist with no mutual interference or is there a compe- 
tition between superconducting and magnetic order? 

Opposite answers come from different experiments. 
Recent muon spin rotation experiments on Ru-1212 2 
have suggested that the magnetic moments are not af- 
fected by the appearance of superconductivity below 
45 K. This assumption motivated Shimahara and Hata 
to study the superconducting properties of a system con- 
sisting of alternating SC and FM layers in the presence 
of a fixed internal magnetic field, neglecting the effect of 



SC on FM&. On the other hand, experiments on chem- 
ical substitution (doping) of Ru-1212 indicate that the 
magnetic and superconducting critical temperatures, T m 
and T c , respectively, are affected in an opposite way with 
decreasing T c and an increase of T m from their refer- 
ence values T c ~ 45 K and T m ~ 132 K in undoped 
Ru-1213M. These results are indicative of a strong com- 
petition between superconductivity and magnetism. 

Since chemical substitution usually affects several pa- 
rameters at the same time, causing changes of the mi- 
crostructure of the sample, a better indication comes 
from the effect of hydrostatic pressure on the super- 
conducting and magnetic phases. New experiments in 
hydrostatic pressure on Ru-12122i u i show that both T c 
and T m increase linearly with pressure P, but at dif- 
ferent rates. This rate is distinctly larger for T m than 
for the superconducting T c , with dT c /dP sw 1 K/GPa 
and dT m /dP sa 6.7 K/GPe 10 . The relatively small pres- 
sure derivative dT c /dP for the superconducting critical 
temperature is interpreted as an immediate consequence 
of a competition of ferromagnetic and superconducting 
phases: a stronger enhancement of the magnetic phase 
results in a reduced pressure effect on T c as compared to 
underdoped high-T c compounds 10 . 

In this paper, we investigate the effect of hydrostatic 
pressure P in a phcnomcnological model of ferromagnetic 
superconductors 11 with two types of carriers pertaining 
to different layers and responsible for superconductivity 
and ferromagnetism, separately. Our main concern is 
the dependence of the critical temperatures T c and T m 
on pressure and filling, that could shed some light on the 
microscopic mechanism of the coexistence or interplay 
of ferromagnetism and superconductivity. The analysis 
follows from a detailed study of the dependence of the 
model parameters (hopping and exchange integrals) on 
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pressure. The paper is organized as follows. In Sec. |n] 
we introduce the microscopic model and discuss the phase 
diagram and the order parameters at finite temperature. 
In Sec. IIIII we give a detailed description of the depen- 
dence of the model parameters on pressure. In Sec. IIVI 
we report the numerical results for T c and T m . Finally, 
in Sec. we present our conclusions and give directions 
for future work. 



II. TWO-BAND MODEL 

The unit cell of Ru-1212 may be described as a 'dou- 
ble bilayer', each bilayer being composed by a Cu02 
and a RuC>2 layer, separated by an intermediate SrO 
layer. The two bilayer blocks are in turn separated by 
a Gd ion, which also serves as an inversion point for 
the unit cell (see, e.g., Fig. 1 in Ref. IT2T ). As in the 
high-T c cuprates, superconductivity is believed to set in 
within the CUO2 layers, while ferromagnetism may be 
thought as mainly due to the ordering of the Ru mo- 
ments in the RuC>2 layers. This has suggested that both 
the SC and FM phases in the ruthenocuprates are not 
homogeneous at the microscopic scale. In particular, the 
SC order parameter may develop^ a spatial variation 
with non-zero total momentum as in the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO or LOFF) phaso 13 i 14 , while 
recent experimental results indicate that the FM order is 
characterized by predominant AFM domains separated 
by nanoscale FM domains^ -. A separate origin of SC 
and FM correlations is also supported by band structure 
calculationsi2ii£, which clearly indicate the existence of 
electronic subbands pertaining to the CuC>2 and to the 
Ru02 layers, respectively, as well as of a hybridization 
term, due to the bridging apical oxygens between adja- 
cent layers. Electronic subbands in the Cu02 and R11O2 
layers are mainly characterized by the Cu-3d x 2_ y 2 and 
Kn-4:d X y orbitals, respectively, as well as by the 0-2p xy 
orbitals±2J£. 

In order to study the coexistence of SC and FM in fer- 
romagnetic metals, such as UGe2 (Ref. yj), a single-band 
model has been originally developed within the mean- 
field approximation by Karchev et ali&, and then nu- 
merically discussed by Jackiewicz et alH. It has been 
pointed out, however, that a single band model does not 
produce coexistence, but rather a first order transition 
between phases^. On the other hand, coexistence of FM 
and SC is permitted in multiband models 19,20 . 

In the case of the ruthenocuprates, a minimal model 
for coexisting FM and SC is then the two-band model of 
Cuoco et al, 11 . There, one may additionally allow for the 
hybridization of the two bands pertaining to the Cu02 
and RuC>2 layers, respectively, by explicitly including an 
interlayer hopping term. Accordingly, we assume that 
the total Hamiltonian can be decomposed as 



Htot = Hsc + Hfm + ffsc-FM, (1) 
with Hqc and Hfm describing the Cu02 and Ru02 lay- 



ers, respectively, while -ffsc-FM contains both the hy- 
bridization term and exchange correlations between the 
two subbands. Within the mean-field approximation, the 
three terms read in turn: 

H sc = ^C k dL^-E[ A 4 T dt -ki+H.c.](2a) 

kcr k 

E(^+r M ) c L<w, (2b) 

kcr ^ ' 

[*±(4 CT Ckcr + H.C.) 
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krr 
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Here, c^, d\ a [ckcr, d\ s _ a \ are creation [annihilation] oper- 
ators for electrons, with wavenumber k and spin projec- 
tion f7 = ±or(7G{|,J,} along a specified direction, in 
the Ru02 (FM) and CuC-2 (SC) subbands, respectively; 
£k = e| c — \i and = e™ — fi are the subbands dis- 
persion relations, measured with respect to the common 
chemical potential fi; t± is the (momentum conserving) 
interlayer hopping or hybridization term between the SC 
and FM subbands; j± = J±/J\\ denotes the ratio of the 
interlayer to the inplane exchange couplings. In Eqs. J5J, 



A = .^(rf-kiit)' 

k 

M = J w H (( c kT c kT> - (cL ki C-ki; 



(3a) 
(3b) 



are the mean-field SC and FM order parameters, respec- 
tively, (. . .) denotes a self-consistent statistical average, 
and g > is the SC coupling constant, which we assume 
to be independent of momentum, for the sake of simplic- 
ity. 

Neutron scattering experiments on the magnetic sus- 
ceptibility of the ruthenocuprates^ agree fairly well with 
a MF picture, thus indicating that a MF desription is ad- 
equate to describe ferromagnetism in the Ru02 planes. 
In order to improve the MF approach, one should take 
into account for the effect of spin density fluctuations via 
a dynamical susceptibility or vertex corrections. This im- 
provement has not been considered in our model, where 
the Stoner criterion has been used in order to describe 
the essential aspects of magnetism. 

For the band dispersions, within the rigid tight-binding 
approximation, we take 



SC,FM 



-2i(cos k x + cos k y ) + At' cos k x cos k y , (4) 



where t, t' are the appropriate nearest neighbor (NN) and 
next-nearest neighbor (NNN) hopping amplitudes for the 
two layers, respectively. 

Apart from a constant term, Eq. Q can be conve- 
niently rewritten in matrix form as 



Htot = Yj B^HB^ 



(5) 
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where — (d^ d_kj, c-ki) is a four-component 
spinor accounting for the two different orderings, the real 
symmetric matrix 



H = 
fa 



-A 

tx 




-a 



-A tx 

f- \]xM -t ± 

Ck + \M 

-*± -Ck + \M) 



(6) 



has been introduced. It is worth noting that H may be 
thought as being composed of four 2x2 blocks, each diag- 
onal block pertaining to the SC and FM subsystems, re- 
spectively. Competition between SC and FM is provided 
not only by the off-diagonal blocks, which only contain 
the interlayer hopping term tx, of kinetic origin, but also 
by the magnetization-induced splitting of the SC sub- 
band, induced by the interlayer exchange coupling J±. 
Inversely, the presence of superconducting correlations 
in the CuC"2 layers (A ^ 0) does not explicitly enter the 
(2, 2) FM block, if not, e.g., through the common chem- 
ical potential fj,, to be self-consistently determined as a 
function of the total number of electrons. 

Eq. |JSJ can be diagonalized by means of standard 
techniques in terms of the four real eigenvalues E^ a 
(a = 1, .. .4) of the band matrix, Eq. JBJ. The SC and 
FM order parameters, A and M, can then be derived 
self-consistently from Eqs. as2i: 



m = 



g \ ^ dE ka ( (3E ka 
OA m I 2 

k,a v 



4N 



J\\ \ - dEk a 

27v -Lw tanl \ 



(7a) 
(7b) 



at fixed number of electrons N and inverse temperature 
(3 = 1/k-oT. Numerical analysis of Eqs. Q shows 11 
that the two-band model allows for the coexistence of 
SC and FM over a reasonable range of parameters. In 
Eqs. 0, the eigenenergies E ka are implicit functions 
of the order parameters A, M, via the secular equation 
det(il — E^al) = 0, and their derivatives can be calcu- 
lated by means of the implicit function theorem (Dini's 
theorem). By direct inspection of the secular equation, 
it can be shown that the eigenvalues E-^ a are even func- 
tions of A, while there always exist paired branches a, a 
such that Eta (~ M) = —E^aiM). The critical temper- 
atures T c and T m for the onset of SC and FM are de- 
fined as the largest temperatures for which Eq. (|7a|l and 
(|7bl) have nonzero solutions A and M, respectively. They 
have been obtained by simultaneously solving Eqs. 0, 
linearized with respect to the appropriate order parame- 
ter. 



III. PRESSURE DEPENDENCE OF THE 
MODEL PARAMETERS 

In Eqs. for the order parameters, pressure P en- 
ters through the band (t, t' in each layers, and tx) 
and the coupling parameters (<?, J», Jx), as well as 
through the doping level, here parametrized by the chem- 
ical potential fi. The phase diagram of correlated sys- 
tems close to an ordering instability is usually character- 
ized by the interplay of a pressure-induced doping vari- 
ation and any other 'intrinsic' pressure effect, here ac- 
counted for by the pressure-dependence of all other model 
parameters^. This scenario can qualitatively explain the 
pressure dependence of T c in the high-T c cupratesS, in 
particular also when an anisotropic doping redistribution 
takes place among inequivalent layers due to an applied 
pressure^. This scenario has been recently related to the 
proximity of an electronic topological transition, where a 
pressure- or strain-induced change of the topology of the 
Fermi surface takes place either because of a change of 
the electronic structure (at constant doping), or because 
of a doping variation (at fixed or rigid band structure, as 
is usually assumed)^. 

An estimate of the pressure dependence of the band 
parameters could in principle be achieved through ex- 
periments or by extensive ab initio calculations 26 . How- 
ever, due to the limited number of experimental results 
on ruthenocuprates in hydrostatic pressure, in the fol- 
lowing we will discuss a simplified scheme allowing us 
to describe the pressure variation of the relevant model 
parameters. 

We will be mainly concerned with the pressure de- 
pendence of the band parameters and of the ex- 
change integrals. Although a pressure dependence of 
the superconducting coupling parameter g is also to 
be expected on general grounds, its actual functional 
form would depend on the microscopic mechanism of 
superconductivity^ 2 ^ 3 ., which is currently a matter of de- 
bate for the ruthenocuprates. In view of the reduced 
pressure effect on T c , as compared to T m , we will then ne- 
glect altogether the pressure dependence of g, although, 
as mentioned in Sec. [fl the enhancement of the magnetic 
phase with increasing pressure could justify a reduction 
oidT c /dP. 



A. Band parameters 

The relevance of the tight-binding approximation for 
modelling the band structure of the Cu02 and the 
Ru02 layers in both cuprate and ruthenate compounds 
has been reviewed by Mishonov and Pene^^. A 
pressure-induced variation of the band parameters en- 
tering Eq. (0J, namely the nearest-neighbor (NN) and 
next-nearest-neighbor (NNN) hopping amplitudes t and 
t', respectively, and of the interlayer hopping amplitude, 
tx, entering Eq. JSJ, can be approximately accounted for 
within the extended Huckel theory^. In this context, 
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such parameters can be roughly approximated by the 
overlap integrals between the appropriate orbitals, which 
are the Cn-3d x 2 _ y 2 and the 0-2p xy , for t sc ; the Ku-4:d xy 
and the 0-2p x , y , for t FM ; the 0(lj-2p x and 0(2)-2p J/ , for 
t' in both layers; and the Ru-4<i3 Z 2_ r 2 and 0-2p z , for t±. 
These are two-center integrals, which have been evalu- 
ated analytically in terms of the distance (and relative 
orientation) of the two orbital centers. At large inter- 
site distances, the approximate behavior employed e.g. 
in Ref. |28| is recovered. 



B. Exchange integrals 

In order to calculate the dependence on hydrostatic 
pressure of the exchange interaction in the RuC>2 planes, 
Jy, and the interlayer exchange coupling Jj_, we fol- 
low the approach of Munro2& using the same approxi- 
mation scheme and generalize it to the case of Ra-4d xy 
and Cu-3rf 2 .2_ l) 2 orbitals. In this approach the quantity 
J~ 1 (dJ/dP) is determined within the theory of solids un- 
der hydrostatic pressure in which the application of pres- 
sure is represented in terms of the crystal compressibility 
and two other parameters associated with the electronic 
screening (A) and the wave-function-distortion (O)iSS To 
a certain extent, therefore, these two parameters take 
into account for the many-body effects in an "equivalent 
pressure-free" model system. 

The generic exchange integral at zero pressure is de- 
fined by: 



dVr^tiYiWji^)— &(r 2 )^(n), (8) 



where r\ 2 = |ri — r2 1 , and ipi(r) — i/j(r — r^) and 
4>j(v) = <p{v — Yj) are the appropriate hydrogenoid or- 
bitals on atoms i and j, respectively. We seek for the 
pressure dependence of the expression (JSJ) which is to 
be approximately determined as a function of A, f2 and 
the compressibility k. To this aim, one first postu- 
lates a scaling of the charge-related coupling constants 
which can be written in the forms e 2 — > A(P)e 2 and 
Z — > il(P)Z/A(P), where Z is the effective charge num- 
ber of the nuclear unit. Second, one assumes that the 
fractional variation of the one-electron state tpi and 4>i 
can be written as a function of k, A, and 17, i.e. 



1 dip 
~^~dP 

6dP 



/(K,A,n), 

g(K,A,fl). 



(9a) 
(9b) 



we can write 



1 d J 1 „ ,./ . „ / . dA 

JdP ~ 3 K+ ' )+ ' )+ dP 



d^riri) / d 3 r 2 ^(v 2 



K 1 

_ 3Jj J r 12 

x [ Tj ■ Vf(r 2 - Ti ) + Tj ■ V0(n - r,)} . (11) 

where we have taken = and the last term comes 
from an expansion around = 0. In evaluating I jllj l. 
we make use of the following approximations 2 ^. First, it 
is assumed that the major contribution to the wavefunc- 
tion distortion comes from its radial part, R(r) say, so 
that ip~ 1 dip/dP w R~ 1 dR/dP. Second, from the same 
assumption it follows that 



' dr 



OR 



(12) 



where p c = Zc/2ao, c is the distance between the ions 
and ao is the Bohr radius. 

In our specific case, to evaluate J» we must consider 
the Ru Ad xy orbitals, whose radial function, assuming a 
hydrogenoid-like wavefunction, is: 



Ridir) 



1 Q7/2 ( JL 
24v/l0 V 2fl o 



3/2 



e - np/2 p 2 {6-n P ), (13) 



where p = Zr/2ao. One thus finds 
1 8R Ad 1 (1 „ p 



Rid dP 



1/7 

2 \n 



Q-VLp 



p 



dn 
dp' 



(14) 



In this expression, one makes use of the approximation 
p pa (p 2 )id/ {p}id instead of the approximation p ps {p)id 
to retain a better numerical accuracy. Using Eq. I|12fl , we 
also have: 

rj .v*-( S (&) -i^-e(^)J** 

(15) 

Setting f — g into Eq. Ijllll and using the relations l|14l - 
I15p . we obtain the following expression for the variation 
of the exchange integral Jii on pressure: 



1 dJ\\ 

~j\\~dT 



K 1 + p c - 6 



Pc 

P I nl 



12 



Pc 



dA 
dP 



2 7 



(p 2 )m 
(P) 



-id 



6-p 



id 



p( 6 - P) / id 
dQ. 
dP 



(16) 



Using the fact that 



where we have expanded up to the lowest contribution 
in Q in the coefficients. This expression is equivalent to 
Eq. (10) of MunroS for the 3d orbitals. 

In the evaluation of Jj_ we must evaluate the exchange 
integral between the Cu 3d x 2_ y 2 orbital (ip) and the Ru 
Ad xy orbital (</>). The radial function for the 3d x 2_ y 2 
orbital is: 



d(l/r 12 ) 
dP 



(10) 



1 ' 9\/30 \2ao 



3/2 
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so that we derive: 

1 8R 3d 1/7 ( P 2 ) 3d \ dn 



R 3d dp 2 \n ( P ) 3d J dp' ' 

Making use of (|12fl and i|17|) . we obtain: 



v r V4>^[2{^-^--p c njR 3dl (19) 

where p c = Zc/2ao, and c is the distance between the 
Ru-Cu ions, this time. The functions / and g are now: 
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2 \0 \6- ftp, 4 „, 
1/7 (p 2 ) M \<m 



(P) 



4(/ 



2\n ( P ) 3d J dp' 



(20) 



Using Eqs. JTUJ) and in Eq. f|XTf> . we obtain the 
following variation of J± with pressure to the lowest order 
in f2: 



1 dJ 



_L 



^ 1 

./ dP~ ~ 3 

7- 



« 1 + p c - 4 



P / ::,/ 



4rf 



<p 2 ; 



.!(/ 



4J 



<P> 



-4 



:!(/ 



6-p 



dA 
dP 

id 



dn 
dp' 



(21) 



As a first step, we assume that the contributions from 
dSl/dP and dA/dP in Eqs. O and are small com- 
pared to the one from the compressibility, and neglect 
them altogether. In so doing, the quantities to determine 
in order to have the full dependence of the exchange in- 
tegrals on pressure, are the compressibility k, p c , p c and 
all the average values that appear in Eqs. (fTf>)) and (|2T|l . 

In determining p c , we need to know the Ru-Ru ions 
distance c/ao- The crystal structure analysis of Ru- 
1213i£ gives for the distance between Ru ions in the RuO 
layer and the apical oxygen of RuC>6 octahedra the value 
d(Ru — O ap i) = 1.912 A. We can then determine the Ru- 
Ru distance in the layer as c = 4d(Ru — O ap i) tan ^ = 
4.415 A, or c/ao — 7.56144. Consequently, we ob- 
tain p c ~ 15.123. From the crystal structure analysis 
we also know the distance between Ru and Cu ions, 
ef(Ru- Cu) = 4.102 A and so we obtain p c ~ 15.508. 
All the average values that appear in Eqs. and ij?T|) 
have been evaluated analytically by the use of the ra- 
dial functions Eqs. i|13|) an d (|17|) . and their values are 
reported in Table [I] 

Finally, we need the compressibility k. To our knowl- 
edge, no experiment has been yet performed to determine 
n. We can estimate this quantity based on the pres- 

9.10 



sure dependence of T m known from experiments****, viz. 
dT m /dP = 6.7 K/GPa. Within Stoner's model of ferro- 
magnetism, one has k^,T m — 7r~ 1 [6(a— l)/ai?] 1 ^ 2 , where 
R = (p'/p) 2 — p" I Pi P — p(p) is the density of states 
(DOS) at the chemical potential p, and a = Jup, a > 1 
being the Stoner criterion for ferromagnetism 3 ^. Neglect- 
ing the pressure dependence of the DOS, one roughly 



id 



(P 2 ) 



1 

Pi h 
P 



6-P/ 4d 

1 



: 10.5 

125.8 
= 0.125 

= -1.75 

= -0.125 



-Pi 4d 
1 

p(6 - P) 



= 



(p)-id 

3d 



(P 2 > 
1 

P 



= 7 
= 56 

0.166 



TABLE I: Average values entering Eq. (1161 and Eq. I21L as 
analytically evaluated by means of Eq. 1131 and Eq. HI 71 1 for 
the radial functions. 



finds 



1 dT m 
T m dP 



1_1 1 gjjj 

2a- ll] l ~dP 



(22) 



which explicitly depends on filling through a. This has 
to be contrasted with the relation^ 



1 dT m 
T m dP 



1 dJ\\ 
~dP' 



(23) 



holding within the Heisenberg model of ferromagnetism 
at the mean-field level. Making use of the latter, albeit 
filling- independent, expression and of Eq. (|16f) . we can es- 
timate the compressibility as k — 3.4 -10 -3 GPa -1 , which 
is a reasonable value, if compared to the values known for 
other cuprate materials with perovskite structur e) 31 ! 32 . 
Inserting this value back in Eqs. (|16|) and 1|21|) . we ob- 
tain the dependence of the exchange integrals on pres- 
sure. Comparing the relative pressure coefficients, we 
also obtain (dlnJii/dP) / (d In J± /dP) ~ 0.77, indicat- 
ing that J± increases faster with pressure as compared 
to J||. We would like to stress that the relation of T m 
with the relevant model parameters could be different if 
approximations beyond the MF level are taken into ac- 
count. Neverthless, ac-susceptibility and resistivity mea- 
surements on ruthenocuprate a 9 i 10 fairly well agree with 
the MF picture, thus indicating that a MF desription is 
adequate to describe ferromagnetism in the Ru02 planes. 



IV. NUMERICAL RESULTS 

In Fig. ^ we show our numerical results for both 
the superconducting critical temperatures, T c , and the 
ferromagnetic critical temperature, T m , as a function 
of the chemical potential p, for 11 values of pressure 
P = — 2 GPa. Fig. [21 shows the shapes of the Fermi sur- 
face relative to the Ru02 layer, e™ = p*, corresponding 
to the chemical potential p* which maximizes Tm at a 
given pressure. One immediately concludes that, within 
the present approximation, pressure has a negligible ef- 
fect on the optimal filling for T m . 

At zero pressure, we take t — 0.3 eV and t'/t — 0.45 for 
both the SC and the FM banda 12 i 15 i 33 , and t± = QML 
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The values of the coupling parameters at zero pressure 
have been chosen so to reproduce the observed optimal 
values of T c and T m at P = 0. Specifically, we take 
g = 0.042 eV, Jy = I At, J± = O.lt. 

At zero pressure, we find three separate regions of co- 
existence of superconductivity and ferromagnetism, with 
T m displaying three pronounced 'domes' as a function 
of chemical potential \i. This is mainly a consequence 
of the Stoner criterion, which in the simplest version, i.e. 
neglecting interlayer exchange and in the absence of com- 
peting SC order, reads J\\p > 1, so that ferromagnetism is 
enhanced where the DOS is largest, i.e. close to Van Hove 
singularities or electronic topological transitions^. As 
pressure increases, the band widens and the DOS peaks 
lowest. However, the exchange couplings J\\ and J± are 
expected to increase, as a result of a larger overlap of the 
orbitals in Eq. (JSJ) , so that the Stoner criterion is satisfied 
over larger filling ranges. As a consequence, separate re- 
gions of coexistence of SC and FM are expected to merge, 
as shown in Fig. ^ As is also shown in Fig. ^ the fer- 
romagnetic transition temperature is found to increase 
with pressure at a rather large rate, in good qualitative 
agreement with experiments. Its values go from 130 K 
to 160 K in the pressure range P = — 2 GPa. Thus 
the ferromagnetic state appears to be strongly stabilized 
under pressure which should have some consequences for 
the superconducting state. 

As mentioned in Sec. lIIII up to now we have neglected 
any explicit pressure effect on the SC coupling parame- 
ter g. This is motivated by a lack of either theoretical or 
phenomenological input for the microscopic mechanism 
of superconductivity in the ruthenocupratcs, which is ex- 
pected to be of unconventional nature, as is possibly for 
the high-T c cuprates. Quantitatively, this approximation 
should not affect much our results, in view of the small 
pressure effect on T c , as compared to T m . However, we 
have numerically studied the competition of SC and FM, 
both at zero pressure and for increasing P, by tentatively 
assuming a small linear dependence of g on pressure. In- 
deed, the effect of a competing FM phase at P — does 
decrease T c , compared to the case Jn = J± = case, as 
already observed by Cuoco et ali^. This tendency is also 
confirmed at nonzero pressure. 

Since we neglected any explicit dependence of the SC 
coupling constant g on pressure, the albeit small increase 
of T c shown in the inset of Fig. ^ must be mainly at- 
tributed to the pressure-induced changes of the kinetic 
terms in the SC and SC+FM Hamiltonians, Eqs. l(2*ajl 
and (|2c|) . i.e. changes in the band structure. Although 
the fine details of the variations of T c and T m are related 
to each other in an inherently nonlinear way through 
Eqs. 10, one expects that a pressure- induced enhance- 
ment of the hopping parameters i, t', and t± results in 
a shift towards the band bottom of the Van Hove sin- 
gularity pertaining to the SC subband, accompanied by 
a steepening of the DOS, which is indeed recovered in 
the tendency of the maxima in the T c curves to move to- 
wards lower chemical potentials with increasing pressure 
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FIG. 1: Superconducting critical temperature T c (solid lines, 
and inset) and ferromagnetic critical temperature T m (dashed 
lines), as a function of chemical potential fi, for different pres- 
sures P = — 2 GPa. Lower curves correspond to lower pres- 
sures, as indicated by the arrow. 

3 p ' 1 1 1 pLo ' : 

P = 2 GPa 
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k x 

FIG. 2: Fermi lines of the RUO2 layers corresponding to the 
maximum T m in Fig. □ for P = and P = 2 GPa. 



(Fig.n inset). 



V. CONCLUSIONS 

We have considered the two-band model for the coex- 
istence of superconductivity and ferromagnetism in the 
ruthenocuprates. We have self-consistently solved the 
equations for the SC (respectively, FM) critical tempera- 
ture in the presence of FM (respectively, SC) order, both 
as a function of filling (here parametrized by the chemical 
potential /x), and as a function of pressure. We find sep- 
arate filling ranges where the coexistence of SC and FM 
is allowed, merging into larger ranges, as the Stoner cri- 
terion gets more effective with increasing pressure. The 
ferromagnetic transition temperature is found to increase 
with pressure at a rate distinctly larger than that of the 
superconducting temperature, in good qualitative agree- 
ment with recent experiments in the ruthenocuprates. 
Due to the competition between superconductivity and 
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fcrromagnctism, the stronger enhancement of the mag- 
netic phase results in a suppression of the pressure effect 
onT c . 

Our present model does not account for spatial varia- 
tions of both order parameters, or of their symmetry in 
k-space, which has been assumed to be s-wave, for the 
sake of simplicity. In particular, a c-axis modulation of 
the order parameters may be important in view of the 
stronger pressure dependence of the interlayer correla- 
tions. A more detailed study of the pressure dependence 
of the critical temperatures would require more reliable 
estimates of the compressibility and of the pressure de- 



pendence of the charge filling, which await more experi- 
mental work. 
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